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The bonding and stability of MgSi and AlMgSi compounds relevant to AlMgSi alloys is investi- 
gated with the use of (L)APW+(lo) DFT calculations. We show that the f3 and j3" phases found 
in the precipitation sequence are characterised by the presence of covalent bonds between Si-Si 
nearest neighbour pairs and covalent /ionic bonds between Mg-Si nearest neighbour pairs. We then 
investigate the stability of two recently discovered precipitate phases, Ul and U2, both containing 
Al in addition to Mg and Si. We show that both phases are characterised by tightly bound Al-Si 
networks, made possible by a transfer of charge from the Mg atoms. 
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I. INTRODUCTION 

Precipitation or age hardened alloys are today one 
of the most important alloy types in industry. In the 
AlMgSi alloy system, Mg-Si and Al-Mg-Si precipitates 
formed during specific heat treatments give rise to a very 
significant increase in strength. The precipitation se- 
quence is generally accepted to be: 

SSSS -> Mg/Si Clusters GPZ -» 0" 0' -+ 0, 

(1) 

where GPZ refers to a Guinier Preston Zone and SSSS 
refers to a Super Saturated Solid Solution. Very little is 
known about the early stages of the precipitation pro- 
cess. However, it is believed that when the SSSS is 
heat treated Mg and Si atoms quickly diffuse substitu- 
tionally to form small clusers due to the large amount 
of quenched-in vacancies^ (a large part of the vacancies 
move to interfaces like surfaces and grain boundaries in 
the later stages of the heat treatment). Although the 
details are hard to investigate experimentally, several 
studies of such clustering have been carried out using 
atom probe microscopy^. The first phase which can 
be resolved using high resolution electron microscopy 
(HREM) is the GPZ. From this, a model of its crystal 
structure has recently been proposed^. The structure of 
the proceeding phase, 0", was also solved using electron 
microscopy techniques^, a result which has later been 
supported by ab initio calculations 6 . 

Earlier it was believed that the GPZ, 0" and phases 
all had the stoichiometry of the phase, Mg2Si, and that 
alloys should be optimized accordingly. Thus, the termi- 
nal phase was of primary importance. Later however, 
it had been confirmed that the 0" gives a greater contri- 
bution to the hardness due to its semi-coherent interface 
with the aluminium matrix and needle-like shape, which 
is more effective for dislocation pinning 7 . The 0" stoi- 
chiometry has been shown to be Mg 5 Sig£. 

Recently, several additional phases have been identi- 



fied experimental!} 
tion sequence: 



giving the extended precipita- 



SSSS -► Mg/Si Clusters -» GPZ 
-> (0' + U2 + U1 + B') 



0" 
>0. 



(2) 



In ref. 0the phases Ul, U2 and B' are referred to as type 
A, B and C respectively. These three phases, in addition 
to 0' , are often grouped together since little is known 
about their interdependence. However it is believed that 
the peak of concentration with respect to time for each of 
these structures follows the ordering given by the above 
precipitation sequence^. They all form relatively late in 
the precipitation sequence usually at temperatures in the 
range 200-300 °C, and in Si-rich alloy compositions. 

It has been considered a general rule of thumb that 
successful aluminium precipitation hardening alloys con- 
tain secondary and ternary alloying elements which are 
larger and smaller than aluminium 11 . Now, the concept 
of atom size in this context must be based on the type of 
bonding involved. One can have either ionic, metallic or 
covalent radii for the constituent elements giving dramat- 
ically different values for atomic sizei^. It is clear that 
when studying the electronic density of compound struc- 
tures this type of concept may lead to an oversimplifica- 
tion. It it therefore interesting to carry out a theoretical 
study of the electronic structure and bonding character- 
istics of the relevant phases with respect to their relative 
stability. This is the purpose of the present work, where 
we employ full potential ab initio methods based on Den- 
sity Functional Theory (DFT) to investigate the bonding 
within the above mentioned structures. In section II wc 
describe the Linear Augmented Plane Wave + local or- 
bitals approach, (L)APW+(lo), used for all calculations. 
Section III describes the results obtained for the three 
models of precipitate phases containing only Mg and Si 
(0", 0', and 0), and section IV deals with the phases 
containing Al, Mg and Si (Ul and U2). 
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II. METHOD 

The ab initio calculations were performed using 
WIEN2k, a program package implementing the full po- 
tential (L)APW+(lo) DFT method^. The Augmented 
Plane Wave approach (APW)+(lo) method 1 ^ differs from 
the LAPW method in the linearization of the APW's. In 
slater's original APW method^ the unit cell is parti- 
tioned into non-overlapping atomic spheres and an inter- 
stitial region. The basis functions are plane waves for the 
interstitial region and radial wave functions within the 
atomic spheres. In the LAPW method the basis func- 
tions inside the spheres are linearized with respect to the 
energy E t : 

^k n (r) = y^X A im,k n ui(r, Ei) + Bi m ^ n ui{r,Ei))Yi m {f) 

lm 

(3) 

where ui{r, E{) are radial functions and Yi m (r) are spher- 
ical harmonics. ui(r, Ei) is the energy derivative of 
ui(r,Ei). Ai m k n and £?/ m fc m are determined by matching 
the above basis to the value and derivative of the plane 
waves for each k-vector at the sphere boundary. An al- 
ternative way for doing the linearization is the APW+lo 
method. One starts with the original APW's and adds 
local orbitals to obtain the variational flexibility in the 
basis functions. 

^fe„0) = ^2(Ai m , kn ui(r, Ei) + (j)i )Y hn (r) (4) 
lm 

4>io = B[ m ui(r, Ei) + Ci m ui(r, Ei) (5) 

At first sight this looks very similar to the LAPW basis. 
However, the B\ m and Ci m are no longer dependent on 
the wave-vector and are determined by the requirement 
that the local orbital is zero at the sphere boundary and 
normalized. The great advantage of this scheme is that 
the calculations converge to results almost identical with 
those of the LAPW method but for dimensioning param- 
eters which effectively leads to a smaller basis se1*i&. The 
WIEN2k code uses a mixed APW+lo/LAPW basis set, 
exploiting the advantages of both methods. 

To make the results for the different structures com- 
parable, we used the same set of APW+(lo) parameters 
for all calculations: R m t = 2.1 Bohr, R m tK ma x — 7 and 
Gmax — 14 Ry 1 / 2 . Here, R m t is the muffin tin radius, 
Kmax is the plane wave cut-off and G ma x is the maxi- 
mum Fourier component of the electron density. For all 
calculations we used the modified tetrahedron method 1 ^ 
for Brillouin zone integrations. All k-point meshes were 
checked for convergence. Thus in general, the highly 
symmetric structures (like bulk Al, Si and Mg) with few 
symmetrically incquivalcnt atoms in the unit cell, require 
a denser k-point mesh than the precipitate phases with 
larger unit cells and a lower symmetry. For the exchange- 
correlation potential we used the Generalized Gradient 
Approximation (GGA) of Perdew et. ali& 




FIG. 1: Conventional unit cell for fluorite Mg2Si, equivalent 
to the (3 phase in the AlMgSi precipitaion sequence. The large 
spheres are Mg atoms and the smaller ones are Si. 



III. MgSi PHASES 
A. (3 

The j3 phase (fluorite Mg2Si) is the terminal equilib- 
rium structure of the precipitation sequence. It has a 
fee primitive unit cell (space group Fm3m (225)), with 
an experimental lattice parameter a=6.39 A 4 . It forms 
precipitates of a plate-like or cubic shape up to 20 /jm 
in diameter. Its interface with the Al matrix is fully 
incoherent. The conventional unit cell, containing 8 Mg 
atoms and 4 Si atoms, is shown in fig-H As can be seen, 
each of the Si atoms in the structure has 8 Mg nearest 
neighbours, giving each Mg atom 4 Si nearest neighbours 
at the same distance. The Si atoms are arranged as an 
fee lattice interpenetrated by a sc Mg lattice. 

Several ab initio studies of the bonding in Fluorite 
Mg2Si have been carried out in the pasl* 1 ^ ^. How- 
ever, to our knowledge, none of these involved the LAPW 
or (L)APW+(lo) DFT method. For our calculations we 
use the experimentally observed lattice constant for this 
phase in the Al matrix: 6.39 A. Performing a volume 
relaxation, this value differed by only 0.25 % from the 
calculated optimized value 6.37 A, and by 0.6 % from the 
experimental value for bulk Mg2Si, 6.338 A22. The calcu- 
lated bulk modulus, derived from a second order Birch fit, 
was 54.3 GPa, compared to 59 GPa from experiment 22 . 
Using the 6.39 A lattice constant the nearest neighbour 
(nn) Mg-Si distance is 2.77 A, the Si-Si mi distance is 
4.51 A and the Mg-Mg nn distance is 3.29 A. 

Past linear-combination-of-atomic-orbitals Hartree 
Fock and DFT pseudopotential calculations^*^! indicate 
considerable charge transfer from the electropostive Mg 
atoms to the electronegative Si atoms resulting in a 
partly ionic Si-Mg bond for the j3 phase. This result 
is also supported by the present work. Fig. [2] displays 
the bonding charge density for the (101) plane of fig. ^ 
We define the bonding charge density as the difference 
between the converged valence charge density from DFT 
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Bonding Charge Density for Fluorite 




2 4 6 8 10 



FIG. 2: Bonding charge density for (101) plane. The two 
central peaks are Mg atoms. The four corner peaks are Si 
atoms. Thick lines represent positive contour lines and thin 
lines represent negative contour lines. 
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FIG. 3: Muffin tin decomposed DOS for the (3 phase. Note 
that the atomic Mg DOS is half the magnitude of the total 
Mg DOS. The vertical line represents the fermi level. 



and the charge density derived from the isolated neutral 
atoms. Thus the bonding charge density indicates the 
charge transfer resulting from the converged electronic 
structure. In fig. the thick lines represent positive 
contours corresponding to charge transfer to the region, 
while the thinner lines represent negative contours 
corresponding to charge transfer away from the region. 
The concentration of thick lines around the Si atoms 
thus displays a net build up of charge around the Si 
atoms counterbalanced by a net reduction of charge 
around the central Mg atoms, indicating that some 
ionicity is indeed at play. 

Fig-EHdisplays the corresponding total and partial elec- 
tronic density of states (DOS) for the (3 phase. Here 
the partial DOS represents the "Muffin tin decomposed" 
DOS for both the Mg and Si atoms, in which the occupied 
states are projected onto the muffin tin eigenvectors (see 
section II) of the particular atom. In this way, some in- 
formation can be gained on the local DOS at each atom, 
however such a procedure is not complete since the in- 
terstitial region cannot be locally resolved. On the other 
hand, the total DOS represents the DOS derived from 
all the "Muffin tin decomposed" DOS and the interstitial 
DOS of the entire computational cell. Note that the cal- 
culation is based on the primitive unit cell, containing 1 
symmetrically inequivalent Si atom and 2 symmetrically 
inequivalent Mg atoms, thus the partial Mg DOS con- 
tains the contribution for 2 Mg atoms, and the Si DOS 
the contribution from 1 Si atom. 

A general feature of the total and partial DOS in fig. EI 
is the two broad bands below the Fermi energy (indi- 
cated by the horizontal line at 0.26 Ry) separated by an 
energy of approximately 0.4 Rydbergs. That the occu- 
pancy is dominated by Si states, where the partial Mg 
DOS represents the two inequivalent Mg atoms, is a fur- 
ther indication of the presence of ionicity whereby the 
Mg has simply donated electrons to the Si. In addition 
this phase has a band-gap. The magnitude of the gap 



at the T-point is 0.13 Ry (1.77 eV), which is somewhat 
smaller than the experimental value of 2.27 eV for a lat- 
tice constant of 6.338 AS 

Figs.^t an d b now display the s, p and d decomposed 
partial DOSs for both Si and the two inequivalent Mg 
atoms. In this case, the total DOS represents the to- 
tal muffin tin decomposed DOS for the particular atom, 
that is, the sum of the s, p and d partial DOSs. Inspec- 
tion of fig. 0Jl reveals that the lower band consists of s 
states where as the upper band consists of predominantly 
p states with a little s state character. For the case of 
the Mg, fig. the DOS below the Fermi level largely 
mimics that of Si apart from both occupied bands being 
of mixed s and p character. 

Insight into the general features of the DOSs shown in 
fig. 0] can be understood from the perspective of begin- 
ning with pure fee Si with a lattice constant of 6.39A, 
giving a rather large nn Si-Si separation of 4.52A. For 
such a system there can be no strong hybridization be- 
tween the atomic valence s and p states, maintaining a 
gap of approximately 0.3 Ry in the corresponding DOS 
(not shown) . The s and p band centers differ by about 
0.4 Ry, which is not so dissimilar from the 0.49 Ry differ- 
ence between the isolated atom 3s and 3p valence states. 
With the addition of the interpenetrating cubic array 
of Mg atoms, both bands broaden reducing the gap to 
about 0.2 Ryd. Thus in fig. the dominant s charac- 
ter of the filled lower band and the dominant p character 
of the upper filled band arises from the strong onsite Si 
orthogonality requirement, whereas in fig. 0] the mixed 
s and p character of the corresponding Mg bands arises 
from the Mg-Si s and p matrix interaction elements. We 
note that the corresponding heights of the Mg total DOS 
are significantly less than the Si DOS, following approx- 
imately the square-root dependence of a metallic DOS. 
Thus, through charge transfer to the Si, the Mg plays the 
role of "strengthening" the Si backbone lattice, provid- 
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FIG. 4: a) Si and b) Mg muffin tin decomposed DOS for the 
/3 phase. Note that for the Si DOS, the s-curve is practically 
superimposed on the totaZ-curve for the lower band in the 
range -0.4 to -0.25 Ry. Scaled square-root curve shown for 
comparison. 



ing an explanation of the origin of the j3 phase band gap 
through the pulling down of the excited s, p and d states 
of atomic Si. 



B. /3" 

The /3"-phase has a base-centred-monoclinic conven- 
tional unit cell (space group C2 m ) with experimental lat- 
tice parameters a — 15.16 A, b = 6.74 A, c = 4.05 A, and 
7ah=105.3° 5 . Fig.[S]displays the conventional cell for two 
viewing orientations. This phase forms precipitates of a 
needle-like shape, typically with a thickness of 30-50 A 
and a length of 300-400 A. The needle length, which is 
along the < 001 > direction of the /3"-phase unit cell, 
runs parallel to the < 001 > direction of the Aluminium 
matrix. For the present work we use the experimentally 
derived lattice parameters and the relaxed atomic posi- 
tions of the inequivalent Si and Mg atoms by minimizing 
the forces by a Newton-Coates procedure with respect 
to the C2/m space group. The resulting inequivalent 
atomic coordinates are listed in table [I] The given space 
group of this phase was first solved by high resolution 




£-0 + % % 



b) 



FIG. 5: Conventional unit cell for (3" seen from a) < 001 > 
direction, and b) slightly away from the < 010 > cartesian 
direction. The large circles represent Mg and the small circles 
Si. 



TABLE I: Relaxed fractional coordinates for the inequivalent 
atomic positions of the f3" phase. 



Atom 


a 


b 


c 


Mgl 


0.0 


0.0 


0.0 


Mg2 


0.346(1) 


0.071(9) 


0.0 


Mg3 


0.421(6) 


0.063(6) 


0.0 


Sil 


0.055(7) 


0.662(7) 


0.0 


Si2 


0.194(4) 


0.250(5) 


0.0 


Si3 


0.209(1) 


0.627(5) 


0.0 



electron microscopy techniques^ and later supported by 
FLAPW DFT calculations 6 .. Fig. © displays the bond- 
ing charge density contour plot for the (002) plane. A 
dominant feature is the concentration of charge between 
the Si nearest neighbours, and to a lesser extent between 
the Si and Mg nearest neighbours, indicating that cova- 
lency is at play in this system^. Such charge transfer 
to the bonding regions originates from the core regions 
of both the Si and Mg atoms, in addition to the homo- 
geneous interstitial region between the Mg atoms. The 
depletion of charge from the Mg is comparable to that 
seen in the /3-phase, indicating that for this system both 
ionicity and covalency is present in the bonding. The 
nearest neighbour distance range in this structure is 2.39- 
2.53 A for Si-Si and 2.61-2.84 A for Si-Mg. This is not 
dissimilar to the nearest neighbour distance of 2.33 A in 
covalently bonded diamond cubic Si and 2.76 A in the 
covalent/ionic equilibrium (3 phase. Moreover, one of the 
inequivalent Si atoms is sp 3 (tetrahedrally) coordinated 
and another is sp 2 coordinated 6 . Fig.[7|displays the Total 
DOS (including the muffin tin and interstitial regions) for 
the (3 phase. Unlike the f3 phase, there is no band gap. 
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Bonding Charge Density for 002 Plane Section 



a) 




FIG. 6: Bonding charge density for (002) plane section of the 
/3" Structure. Thick contour lines represent a positive charge 
transfer, while thin lines represent a negative charge transfer. 
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FIG. 8: Partial DOS for the a) Si2 and b) Mg3 (see table. □ 
atom in the f3" phase. 



TABLE II: Fractional coordinates for the inequivalent atomic 
positions of the (71 phase 9 . 



Atom 


a 


b 


c 


Mg 


0.0 


0.0 


0.0 


Al 


1/3 


2/3 


0.632(8) 


Si 


1/3 


2/3 


0.243(8) 



However, depressions can be identified at -0.1 and 0.4 
Rydbergs reminiscent of the gaps seen in fig. [3J By in- 
spection of the partial DOS for each inequivalent atom 
we find that for the Si atoms there is a dominant s char- 
acter for the lower energy states (< —0.2 Ry), whereas, 
at about 0.2 Ry, the p character dominates. In between 
these regimes, a mixture of p and s character exists, in- 
dicating strong hybridization. Indeed fig. displays the 
partial DOS for Si2, which from past work is locally sp 3 
coordinated 6 , showing approximate correspondence with 
the diamond cubic Si DOS 24 . In particular the dominant 
s character of the L2' , the s and p character of the Li 
and the dominant p character of X4 the high symmetry 
points. Alternatively, fig. [SJj shows the partial DOS for 
Mg3, which from charge density profile analysis 6 is in 
a homogeneous metallic environment, displaying an ap- 
proximate shifted metallic square-root behaviour. 



IV. AlMgSi PHASES 
A. [71 

The (71 conventional/primitive unit cell is trigonal 
(space group Ps m i) with experimentally derived lattice 
parameters a = b = 4.05 A, and c = 6.74 A. It contains 
1 Mg, 2 Al, and 2 Si atoms giving the formula MgAl2Si2- 
The unit cell and atomic positions are shown in table ITT1 
and fig.|U Experiments show that the < 001 > direction 
of the Ul unit cell runs parallel to the < 310 > direc- 
tion of the fee Al matrbiii. The Ul phase usually forms 
rod-like precipitates with a length of 50-500 nm and a 
width of approximately 50 nm. Performing a volume re- 
laxation, we found the calculated unit cell equilibrium 
volume to differ by only 1% from the experimentally de- 
rived unit cell volume. From a second order Birch fit we 
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FIG. 9: Ul conventional unit cell in a) perspective and b) 
< 001 > direction. The small, black spheres are Si atoms, the 
large dark spheres Al atoms, and the large white spheres Mg 
atoms. 




FIG. 10: U 1 bonding network 



Bonding Charge Density for 1 10 Plane Section 
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FIG. 11: Ul bonding charge density for the (110) plane. 
Thick lines represent contour levels with a positive value, thin 
lines a negative value. 



obtained a bulk modulus of 71 GPa. 

The Ul phase can be categorized as belonging to a 
class of structures given the name CaAl2Si2-type Zintl 
compounds 25 . In Zintl compounds, the electropositive 
elements are thought of as merely electron donors to the 
electronegative elements which thereby are able to fulfill 
the octet rule. For MgAl2Si2 this implies that each Mg 
atom donates two electrons to Al2Si2 for each unit cell. 
This gives A^Si^ - , or two units of AISi" which is iso- 
electronic to A1N, fulfilling the octet rule. The resulting 
layered structure is sketched in fig. EH AISi" constitutes 
a double layer of tightly bound "chair-like" six-membered 
rings separated by layers of hexagonal Mg 2+ . In the Al- 
Si network, each Si atom bonds to 4 Al atoms forming 
an umbrella-like structure, while each Al atom forms the 
more common tetrahedral structure with 4 nearest neigh- 
bour Si atoms The Al-Si bond length within each AISi 
layer is 2.48 A, while the length of the Al-Si bonds con- 
necting the two AlSi~ layers is 2.62 A. The length of the 
Mg-Si bond is 2.86 A (see below and fig. ITTj) . somewhat 
larger than the length of the partly ionic bond in the (3 
equilibrium structure, but similar to Mg-Si bonds found 
in the Mg 5 Si 6 f3" phased 

The calculated bonding charge density for the (110) 
plane is shown in fig 1111 There is a clear concentration 
of charge between the Al and Si atoms indicating strong 



Al-Si bonding. There is also a small build-up of charge 
between the Si and corner Mg atoms which constitutes a 
coupling between Mg 2+ and AISi" layers. 

Evidence of the layered structure can also found in the 
DOS. Figs. 1121 and 1131 show the muffin tin decomposed 
total and partial DOS. The total DOS below the fermi 
level is characterized by two broad bands, separated by 
a 0.1 Ry band gap (fig. 112b,). It is also evident that the 
major contribution to the total DOS stems from the Al 
and Si states. This is an indication that the Mg do- 
nates charge to the Al-Si network and is not strongly 
involved in the electronic bonding. Furthermore, the Si 
states give a relatively larger contribution to the total 
DOS than the Al states. From the angular decomposed 
Si DOS (fig. 112b ) one can see that the contribution to 
the lower band comes mainly from Si s states, while the 
higher band consists mainly of p states with a modest 
mixing in of s-states. This can be seen as an indication 
of a partial hybridization of the Si states resulting in the 
Si bonding environment being partly covalent. This hy- 
bridization can also be identified in the Al DOS (fig. 113b). 
Compared to a metallic DOS 12 , there is a peak concen- 
tration of s states with mixed in p states slightly below 
Ry resembling the Li lobe of a covalent DOS. However 
for both Al and Mg (fig. 113b ) , the overall shape of the to- 
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FIG. 12: a) total and b) Si muffin tin decomposed DOS for 
the Ul phase. Note that the Al and Si DOS consists of the 
contribution from 2 symmetrically equivalent atoms, while 
the Mg DOS stems from only one atom. One therefore has to 
divide by a factor two to get the atomic DOS for Al and Si. 
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FIG. 13: a) Al and b) Mg muffin tin decomposed DOS for the 
Ul phase. 



B. U2 



tal atomic DOS seems roughly to follow the square-root 
relationship characteristic of a metallic DOS. 

The origin of the band gap at the centre of the occupied 
states can be understood from the same line of reasoning 
as for the (3 phase. In the Ul phase, the Si atoms are 
arranged as an hep lattice with lattice parameters a = 
4.05 A and c = 6.74 A. If we once again consider only the 
Si system (with a nearest neighbour distance of 4.05 A), 
then little hybridization of the atomic Si s and p valence 
states can be expected and a similar DOS to the (3 is seen. 
That is, the occupied valence states will be separated into 
two relatively narrow bands. By introducing the Al and 
Mg, the Al forms a tightly bound bonding network with 
the silicon via the donation of charge from the Mg atoms. 
This coupling results in a broadening (with respect to 
the artificial Si sub-structure) of the two occupied bands 
and in the removal of the band gap seen in [3. This 
has its origins in the stronger mixed s and p character 
of the second occupied band (centred at w —0.2 Ry) of 
Si ffig.ll2bh when compared to the strong p character of 
the corresponding band for (3 (fig. 0} — indicating strong 
hybridization between the Si p states and the Al s and p 
states. 



The U2 primitive/conventional unit cell is orthorhom- 
bic (space group Prima) with experimentally derived lat- 
tice parameters a = 6.75A, b = 4.05A, and c = 7.94A 
It contains 4 Mg, 4 Al, and 4 Si atoms giving the for- 
mula Mg4Ai4Si4 ffig. 1140. The coordinates of the sym- 
metrically inequivalcnt atoms of the unit cell are shown 
in table II I II The morphology and size of the precipitate 
formed by the U2 phase is usually the same as for the Ul 
phase. The < 100 > direction of the unit cell is oriented 
parallel to the < 310 > direction of the fee Al matrh#i£. 
By performing a volume relaxation of the unit cell, we 
get an equilibrium volume differing by less than 1% from 
the experimentally derived unit cell volume. A second 
order Birch fit gives a bulk modulus of 69.1 GPa, very 
close to the value obtained for the Ul phase. 

The U2 phase is similar to the TiNiSi structure 
typo 26 i 27 . As for the Ul phase, the bonding in the U2 
phase is characterized by the electropositive Mg atoms 
donating charge to the electronegative Al and Si atoms 
forming a tightly bound bonding network ffig. H5|) . How- 
ever, unlike for the Ul phase, the electropostitive Mg 
atoms retain a large amount of charge and (using a simple 
chemical picture) there is not enough transfer for the elec- 
tronegative elements to fulfill the octet rule2£. The Al-Si 
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FIG. 15: U2 bonding network 



bond lengths are all 2.59 to 2.61A, while the bond lengths 
of the Mg-Si neighbour coupling is 2.78 and 2.86A 

Fig. 1161 displays the bonding charge density maps for 
the (040), (020) and (040) planes of the U2 unit cell. 
As for the Ul phase, there is a concentration of charge 
between Al and Si nearest neighbours making up the 
AISi bonding network. This charge concentration ap- 
pears to have more of a covalent, directional character, 
than the Ul Al-Si bond. Furthermore, one can identify a 
more pronounced charge concentration in between Mg-Si 
neighbours than for the Ul phase. This fits well with 



TABLE III: Fractional coordinates for the inequivalent 
atomic positions of the U2 phase unit cell—. 



Atom 


a 


b 


c 


Mg 


0.034(9) 


3/4 


0.327(4) 


Al 


0.361(4) 


1/4 


0.432(5) 


Si 


0.239(3) 


1/4 


0.120(9) 



Bonding Charge Density for 040 Plane Section 




a) 2 4 6 8 10 12 14 
Bonding Charge Density for 020 Plane Section 




b) 2 4 6 8 10 12 14 
Bonding Charge Density for 0-40 Plane Section 




FIG. 16: Bonding charge density for a) (040) plane b) (020) 
plane and c) (040) plane of the U2 unit cell. Thick lines 
represent contour levels with a positive value, thin lines a 
negative value. 

the general assumption for this structure class — that a 
large amount of Mg charge is retained. As for the Ul 
phase there is no indication of charge build-up between 
Al and Mg neigbours, indicating that the AISi network 
is mainly coupled to the Mg atoms via the Si atoms. 

This picture is further justified by the total and par- 
tial DOS shown in figs. 1171 and 1181 Comparing the DOS 
of the U2 phase with the DOS of the Ul phase one can 
identify several common features, supporting the propo- 
sition that the electronic bonding picture is very similar 
for the two phases. Note that the total DOS is the sum 
of 1 Mg, 1 Al, and 1 Si atom, as opposed to 1 Mg, 2 Al, 
and 2 Si atom for the Ul phase ffig. H2b l. From the atom 
decomposed DOS of the U2 phase (fig. \17h ) one can see 
that the Al and Si states are dominating with respect 
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FIG. 17: a) total and b) Si muffin tin decomposed DOS for 
the (72 phase. 
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FIG. 18: a) Al and b) Mg muffin tin decomposed DOS for 
the U2 phase. 



to the Mg states, and that the Si states give a slightly 
larger contribution than the Al states. This is consistent 
with the donation of charge from Mg to Al+Si atoms, as 
in the case of Ul. However, the Mg DOS gives a larger 
contribution to the total DOS than for the Ul phase. 
This is an indication that there is less Mg to Al and Si 
charge transfer, than for the Ul Zintl-type phase. We re- 
mind the reader that it is difficult to quantify the charge 
transfer because the (L)APW+(lo) method does not use 
an atomic/localised basis set. Furthermore, the magni- 
tude of the band gap separating the two occupied bands 
is slightly smaller for the U2 phase than for the U 1 phase 
(0.075 Ry vs. 0.1 Ry for the Ul phase). Consequently, 
the region of hybridization around Ry is increased for 
the Al and Si DOS (fig. ED a an d b) consistent with the 
pronounced covalent character of the Al-Si bonds found 
in the bonding charge density. 

As for the Ul and (3 phase, the origin of the band gap 
in the occupied states can be explained by the separation 
of the Si atoms with respect to each other, which for this 
structure forms a slightly distorted h.c.p structure (not 
shown here) with lattice parameters comparable to those 
of the Ul Si hep sub-lattice. 



V. DISCUSSION AND CONCLUDING 
REMARKS 

We have made a comparative study of the electronic 
structure of the (3, j3", Ul, and U2 phases in the AlMgSi 
alloy precipitation sequence. The bonding in the (3 phase 
is characterized by the covalent bonding between Si-Si 
nearest neighbour pairs and ionic/covalent bonding be- 
tween Si-Mg nearest neighbour pairs. For the (3 phase the 
bonding is dominated by the partly ionic Mg-Si bond. By 
calculating the heat of formation for the various precipi- 
tate phases further insight can be gained in the observed 
precipitation sequence. Presently we employ the formula: 

AH = -EAiMgSi - x M gE M g - x S[ E si - xaiE A \, (6) 

where AH is the formation energy /enthalpy per atom, 
^AlMgSi is the energy of the given AlMgSi compound, 
and i?Mg, Egi, and Eai are the equilibrium ground state 
(zero temperature) energies per atom at of hep Mg, dc 
Si, and fee Al respectively. Xi is the relative content of el- 
ement i in the compound. Results for the various phases 
in the precipitation sequence are shown in table IIVI and 
figure El For the (3 and (3" phases we used the exper- 
imentally observed lattice parameters. For the (3" we 
used the relaxed coordinates found from force minimiza- 
tion (table HJ). We have also included the corresponding 
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FIG. 19: Calculated energies of formation for the various 
precipitate phases. The error bars indicate the estimated level 
of accuracy of ±2.0 mRy. The line connecting j3" and j3 is a 
guide to the eye. 



energy from a past model for the f3' phase proposed by 
Matsuda et al£. 

It is difficult to draw any quantitative conclusions 
based on the relative formation energy of these phases in 
Aluminium, at finite temperatures, without taking into 
account the effects of the interface and entropy. The cal- 
culated energies are nonetheless very reasonable. As can 
be seen, only 3 phases (/3, Ul and U2), give a negative 
formation energy. The Ul and U2 are lower in energy 
than (3", but higher than the f3 phase. f3 shows the most 
negative energy, which is line with expectations, since it 
is the terminal equilibrium structure of the precipitation 
sequence. The energy of j3" is slightly positive, but not 
markedly different from (3, Ul and U2. j3", which can 
be derived from an fee super-cell, forms early in the pre- 
cipitation sequence at temperatures above 150 °C, but 
gradually disappears after further heat treatment^. 

We note the unusually high heat of formation for Mat- 
suda's j3' model which has a hexagonal unit cell with 
space group P62m with experimentally derived lattice 
parameters a — b = 4.05A, c = 6.74A and the stoichiom- 
etry Mg2Si. For the calculated formation energy, we used 
the experimentally derived lattice parameters and pub- 
lished coordinates. Upon volume relaxation, we however 
found that the structure reduced in volume by 26% corre- 
sponding to a heat of formation of 25 mRy. Furthermore, 
from a second order Birch fit, we find that for the experi- 
mental lattice constant, the Matsuda's f3' phase is under 
a negative pressure of approximately 6.9 GPa. Since the 
experimentally derived unit cell parameters are regarded 
as being extremely accurate it becomes difficult to see 



how such as phase might exist within the Al matrix. 

To this date U 1 and U2 are the only phases in the pre- 
cipitation sequence of the AlMgSi alloy system containing 
Aluminium in addition to Magnesium and Silicon. The 
crystallographic structure type of both phases belong 
to structure families containing an abundance of doc- 
umented compounds. A characteristic for both phases 
is the formation of a tightly bound network of Al and 
Si atoms made possible by at charge transfer from the 
electropositive Mg atoms. Considering the wide range 
of structures belonging to this family it is possible that 
these these type of structures containing other harden- 
ing elements than Mg and Si would form stable precip- 
itates in Aluminium alloys. For example, it is quite re- 
markable that the Ul phase with Magnesium replaced by 
the rare earth element Europium produces a CaAl2Si2- 
type structure with lattice parameters and coordinates 
differing from the iJl-phase by only few percent^. Sug- 
gested further work include the calculations on these two 
phases with iso-valent species substituted for the Si and 
Mg sites. In the future these calculations could also be 
used together with interfacial energies and entropy cal- 
culations to estimate the solid solubility of these phases 
in Aluminium. 

Furthermore, the rule of thumb that successful alu- 
minium precipitation hardened alloys have secondary and 
ternary elements that are larger and smaller than alu- 
minium is not sufficiently accurate when applied to the 
stability of bulk phases in the precipitation sequence. 
The atomic radii of Al, Si, and Mg vary with the bond- 
ing environment. For example, in the equilibrium j3 
phase the atomic size of Si is increased and that of 
Mg correspondingly decreased in forming a dominating 
ionic/covalent Mg-Si bond. Also, using the appropriate 
values for ionic radiii^, the Mg-Si bond length in the j3 
and P" phases are predicted to be considerably smaller 
than what the present calculations show. A more precise 
formulation would involve the ability of the secondary(A) 
and ternary(B) elements to form stable bonds (A-B, A-A, 
and B-B) with bond lengths which are comparable that 
of the Al-Al bond. 
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Structure 


Space 
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Bulk Modulus 


Energy / Atom 
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Parameters (A) 
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